Detection of Hot Gas in Galaxy Groups via the Sunyaev-Zel'dovich Effect 
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Motivated by the observed shortfall of baryons in the local universe, we investigate the ability of 
high resolution cosmic microwave background (CMB) experiments to detect hot gas in the outer 
regions of nearby group halos. We construct hot gas models with the gas in hydrostatic equilibrium 
with the dark matter and described by a polytropic equation of state. We also consider models 
that add entropy to the gas in line with constraints from X-ray observations. We calculate the 
thermal Sunyaev-Zel'dovich (SZ) signal in these halos and compare it to the anticipated sensitivities 
of forthcoming SZ survey experiments such as ACT, PLANCK and SPT. Using a multi-frequency 
Wiener filter we derive SZ detectability limits as a function of halo mass and redshift in the presence 
of galactic and extragalactic foregrounds and the CMB. We find that group-sized halos with virial 
masses below 1O 14 M can be detected at z < 0.05 with the threshold mass dropping to 3 — 4x 10 13 M Q 
at z 0.01. The SZ distortion of nearby group-sized halos can thus be mapped out to the virial radius 
by these CMB experiments, beyond the sensitivity limits of X-ray observations. These measurements 
will provide a unique probe of hot gas in the outer regions of group halos, shedding insight into 
the local census of baryons and the injection of entropy into the intragroup medium from non- 
gravitational feedback. 



I. INTRODUCTION 

Through the distinctive signatures that cosmological 
parameters, such as the baryon density, matter den- 
sity and spatial curvature have on the cosmic microwave 
background (CMB) anisotropy, recent measurements of 
the CMB temperature and polarization spectra [l], and 
references therein] have placed precise constraints on 
these parameters, thereby establishing a cosmological 
model that is consistent with a wide range of astronom- 
ical observations @, and references therein]. Moreover, 
because the physics underlying the CMB anisotropy is 
to a good approximation linear on these scales, the CMB 
provides a powerful probe into the physics of the early 
universe and the primordial perturbations 0, and refer- 
ences therein]. 

Whereas on large angular scales the CMB anisotropy 
is fairly well described, apart from a late-time integrated 
Sachs- Wolfe (ISW) contribution by the geometri- 
cal projection of inhomogeneities at the last scattering 
surface onto our celestial sphere, the so-called primary 
anisotropies, on smaller angular scales of around an ar- 
cminute CMB photons are more significantly affected by 
gravitational and scattering processes along the line of 
sight. These secondary anisotropies provide an effective 
probe of the physics of the low-redshift universe. 

The most significant secondary anisotropy on scales of 
around a few arcminutes comes from inverse Compton 
scattering of CMB photons off hot gas in galaxy clusters 
along the line of sight, the so-called thermal Sunyaev- 
Zel'dovich (SZ) effect |5j. The population of hot elec- 
trons in the cluster gas imparts energy to photons in the 
Rayleigh-Jeans part of the spectrum pushing them into 



the Wien tail. This distorts the thermal CMB spectrum 
by creating a deficit of photons at low frequencies and 
an excess of photons at high frequencies with no distor- 
tion at the SZ null frequency, v ~ 218 GHz. The galaxy 
cluster appears as a cold spot on the sky at frequen- 
cies below the SZ null and as a hot spot above the SZ 
null. The unique spectral signature of the SZ effect will 
allow multi-frequency, high-resolution CMB experiments 
to make a relatively clean separation of the thermal SZ 
effect from the primary CMB and other contaminants 
such as point sources and galactic foregrounds. 

The utility of the thermal SZ effect as a cosmological 
probe arises from the fact that the Compton distortion of 
the CMB is a scattering effect, so that the central decre- 
ment towards a cluster is independent of the redshift of 
the cluster. Moreover the angular diameter distance flat- 
tens out at z ~ 1 so that the angular size of the cluster 
is approximately constant at high redshift. This means 
that the SZ effect does not suffer from the strong redshift 
dimming of its optical and X-ray counterparts with the 
consequence that a microwave background SZ cluster sur- 
vey can detect a larger proportion of lower mass clusters 
out to higher redshifts. The thermal SZ effect can thus be 
used to construct galaxy cluster catalogs with a well de- 
fined selection function, thereby providing a potentially 
powerful probe of cosmology through the evolution of the 
cluster abundance (for a review see @). 

While most work has focussed on the detection of the 
SZ effect in galaxy clusters and its cosmological applica- 
tion, there is an intriguing possibility that high resolution 
CMB experiments with sufficient sensitivity may be able 
to detect hot gas in lower mass galaxy or galaxy group 
halos. Whereas the central SZ distortion scales roughly 
in proportion to the mass of the halo, the total SZ flux, 
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5sz oc M 5 / 3 /d\, has an additional dependence on the 
angular diameter distance of the halo. This means that 
nearby low mass halos with a large angular size could 
produce a significant SZ flux. This fact was exploited 
in Q to study the possibility of detecting the SZ effect 
in the halo of our neighbouring galaxy, M31, with the 
upcoming PLANCK surveyor. In this paper we focus 
on how well current and forthcoming CMB experiments 
will detect the SZ effect in nearby galaxy and group ha- 
los. Such a detection will provide unique insights into 
the distribution of baryons and the properties of hot gas 
in the outskirts of galaxy and group halos. 

Indeed, in a census of baryons in the local universe 
d, it has been argued that most of the uncertainty 
in the baryon budget comes from the uncertainty in the 
mass in ionized plasma associated with groups of galax- 
ies. Whereas the baryon fraction inferred from light ele- 
ment abundances (e.g., the CMB anisotropy @ and 
the high redshift Lya forest [ll[ give a consistent baryon 
fraction of as 0.044, only about ten percent of these 
baryons are observed to be in stars, hot gas in clusters, 
and neutral and molecular gas in galaxies at low redshift. 
Cool plasma in and around groups of galaxies indirectly 
detected through quasar absorption lines [HI, E3| makes 
up twenty to twenty five percent of the baryon budget, 
while the warm-hot gas residing in the filamentary large- 
scale structure is believed to account for another thirty 
to thirty five percent of the baryons [3, EH- The remain- 
ing thirty to forty percent of baryons is likely to be tied 
up in galaxy groups in a low density, warm-hot plasma. 
Recent X-ray observations [IB] that report average gas 
fractions, / gas « 0.08, out to r$oo in galaxy groups may 
have detected roughly half of these baryons. When com- 
bined with the fact that the gas fraction profiles inferred 
from these measurements are still rising at rsoo (see also 
[13]) this indicates that a large fraction, perhaps as much 
as fifteen to twenty percent, of baryons in the universe 
may be in the outskirts of galaxy halos and in the in- 
tragroup medium beyond rsoo- Recent evidence for this 
warm-hot component has come from quasar absorption 
lines in the ultraviolet (O VL_e.g., [3 01, [II ) and in X- 
rays (OVII and OVIII; e.g., [21() observed around galax- 
ies and groups of galaxies but these detections only map 
the two-dimensional gas distribution along a few lines of 
sight to allow a complete study of its spatial properties. 

More detailed observations of this warm-hot compo- 
nent will yield important clues into the impact of galac- 
tic feedback on the distribution of gas in the intragroup 
medium. The relative dearth of hot gas observed in the 
central regions of galaxies and galaxy groups results in 
these halos being less X-ray luminous for a given tem- 
perature as compared to more massive galaxy clusters. 
The resulting break in the luminosity-temperature rela- 
tion observed in X-ray clusters and groups [22[ suggests 
that non-gravitational processes that modify the entropy 
structure in the central regions of low mass clusters and 
groups are responsible for departures from the cluster 
scaling relations expected in self-similar models. Several 



models have been proposed to explain the excess entropy 
(for a recent review see [23[), including pre- heating of 
infalling gas before it was shock-heated to the virial tem- 
perature, radiative cooling of low entropy gas to form 
stars, and feedback from supernovae and active galac- 
tic nuclei that increased the temperature and reduced 
the density of gas in the central regions of galaxies and 
groups. While a combination of radiative cooling and 
feedback seems to reproduce the central entropy excess 
in low mass clusters and groups, it is clear that improved 
measurements of the properties of hot gas in these halos 
are necessary to understand the exact details of entropy 
injection. Whereas X-ray observations lack sensitivity to 
the hot gas beyond r 50 o because the X-ray luminosity 
scales as the square of the (decreasing) gas density, the 
SZ effect scales linearly with the density and thus pro- 
vides a potentially more sensitive probe of the hot gas 
in the outskirts of galaxy and group halos. Moreover, 
the SZ effect measures the projected gas pressure which 
is a complementary probe to the X-ray luminosity and 
can provide new constraints on models of the intraclus- 
ter and intragroup medium. While the SZ effect has been 
detected in a number of galaxy clusters by single-dish 
and interferometric experiments (for a review see Q) it 
has not not yet been convincingly detected in any galaxy 
groups. However, a new generation of high resolution 
CMB experiments with superior detector sensitivity of- 
fers the hope of detecting and constraining the properties 
of the hot gas in galaxy groups. 

The outline of this paper is as follows. In Section 
2 we construct models of the hot gas in virialized ha- 
los of a given mass and redshift, over the mass range 
M vir = 10 13 M Q - 1O 15 M and redshift range z = 0- 1.5. 
The models take into account constraints on the tem- 
perature, entropy injection and gas fraction from X-ray 
observations. In Section 3 we utilise multi- frequency fil- 
tering techniques to determine the detectability of hot 
gas in nearby halos with current and forthcoming mi- 
crowave background experiments such as the Atacama 
Cosmology Telescope (ACT 0]), the South Pole Tele- 
scope (SPT (H) and the PLANCK mission 26]. In a 
final section we summarize our findings and discuss how 
detection of the SZ effect in group halos will constrain gas 
models and allow a measurement of physical parameters 
such as the entropy injection and baryon fraction in these 
halos. Throughout this paper we assume a flat ACDM 
cosmological model with parameters h = 0.72, Q, m — 
0.26, n b = 0.044, n dc = 0.74, w dc = -1.0, and a 8 = 0.8 
that provide a best fit to the WMAP 5 year data 0. 



II. HOT GAS HALO MODELS 

Historically the isothermal model [27| has been used 
to describe the spatial distribution of hot gas in galaxy 
clusters, primarily to model the X-ray emission originat- 
ing from thermal bremsstrahlung of the hot intracluster 
gas [28J. The (3 model provides a convenient analyti- 
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cal form that has been popular for X-ray surface bright- 
ness profile fitting. However, over the past decade it has 
been realised that X-ray observations of density profiles 
at large radii do not favour the j3 model: the isothermal 
model provides a poor description of the temperature 
profiles of intracluster gas, with observed temperature 
profiles declining at large radii and cooling flows observed 
in the central regions of some clusters [17]. Furthermore, 
it has been shown that there is a fundamental incompati- 
bility between the (3 model parameters fitted using X-ray 
data, and those fitted using data based on the SZ effect 
[29l ] , indicating that the simple (3 model is not sufficiently 
realistic to describe the observed cluster gas physics. 

In this paper, we will study two different models for 
the hot gas in dark matter halos. Our first model, which 
we will refer to as the polytropic model, assumes that the 
hot gas follows a polytropic equation of state and traces 
the dark matter in the outskirts of the halo (e.g. (30l.[3l|). 
As discussed in Komatsu and Seljak [30j this model is in 
good agreement with observed X-ray surface brightness 
profiles and the mass-temperature relation above temper- 
atures of a few keV. Our second model, which we refer to 
as the entropy model, builds upon the first and attempts 
to account for non-gravitational feedback by adding en- 
tropy to the hot gas, similar to Voit et al. 32| . In both 
cases, the hot gas is in hydrostatic equilibrium with the 
underlying dark matter potential. 



A. Dark matter halo 

We define a dark matter halo at redshift z with virial 
mass M v j r and radius r v j r to have a characteristic average 
density equal to A v ; r (z) times the critical density, 
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The Hubble parameter relative to its present value is 
given by 



H(z) 
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where Sl m and fide represent the density of the matter 
and dark energy components, respectively, relative to the 
critical density today. 

In the spherical collapse model [33] for a SCDM cos- 
mology with f2 m = 1 and Qde = 0, the virial collapse 
factor has a constant value A v ; r (z) = 187r 2 . Using nu- 
merical simulations, Bryan and Norman [34j found for 
ACDM cosmologies that the parametric form 



A vir (z) = 18tt 2 + 82[n(» - 1] - 39[n(» - l] 2 



where 
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provides a good fit over a wide range of ACDM cosmo- 
logical models. Once the characteristic average density 
is chosen, the virial mass and radius become uniquely 
related. 

We assume that the dark matter density follows the 
self-similar NFW profile [3a ]. 

Ps 
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where r s is the scale radius of the halo. By integrating the 
density profile and equating the mass found within the 
virial radius to Af v i r , we obtain the density normalization, 
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47rrA r m(c) ' 
where the function, 

m(x) = ln(l + x)- —J— , (7) 
1 + X 

is the dimensionless dark matter mass profile. The con- 
centration parameter is defined as c = r V i r /r s and we 
adopt the fitting formula, 
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to describe the dependence on virial mass and redshift 
following Dolag et al. [3(| . The dependence on cosmolog- 
ical parameters is captured by scaling the virial mass by 
the nonlinear mass M*o- In the spherical collapse model, 
the nonlinear mass is defined as the mass M = (4/3)7rpi? 3 
enclosed within a sphere of radius R for which the vari- 
ance of the linear density field S, smoothed by a tophat 
filter, equals the square of the critical overdensity thresh- 
old S c = 1.68. For the chosen cosmology, the nonlinear 
mass is M* = 2.818 x 10 12 /i~ 1 M Q at redshift z = 0. 

The exact dependence of the concentration parame- 
ter on mass and redshift is still uncertain as various 
parametrizations have been used in the literature (e.g. 
[35l [37l . HH). However, as previously observed [3l|, we 
find that our results are insensitive to the exact choice 
of the concentration parameter. This is due to the fact 
that the SZ effect is less sensitive to the central regions, 
where the changes in the concentration parameter are 
most important. 

During the formation of virialized halos, the collision- 
less dark matter undergoes violent relaxation and the 
collisional baryons get shockhcated. According to the 
virial theorem, the internal energy of a virialized halo is 
twice its gravitational potential energy. From this rela- 
tion, we can define a characteristic velocity dispersion for 
the dark matter and a virial temperature, 
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for the shockheated gas. We assume a fully ionized gas 
with a mean molecular weight /1 = 3+ § x = 0.588 for a 
hydrogen mass fraction of Xh = 0.76. 
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B. Polytropic model 



We first consider a model where the gas follows a poly- 
tropic equation of state, P g oc pj, with polytropic in- 
dex 7. The gas density and temperature profiles can be 
parametrized as, 



Pg( x ) = Pg,ay P oi y {x) , 



T g (x) = T gfl yl Q y y {x) 



(10) 
(11) 



where y vo \y{x) is the dimensionless gas density profile 
and the coefficients p g .o = p g (0) and T g .o = T g (0) are 
two boundary conditions. Assuming that the gas is in 
hydrostatic equilibrium with an NFW potential, we ob- 
tain the analytical solution (e.g. [30|). 



y P o\ y (x) = [1 - Bf(x] 
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For a given mass and redshift, there are 3 free parame- 
ters: the polytropic index 7, the central gas temperature 
T g fi, and the central gas density p g ,o- I n Komatsu and 
Seljak [3^, Hl[ these parameters were specified by requir- 
ing that the gas density profile matches the dark matter 
density profile in the outer parts of the halo. This as- 
sumption is known to be in good agreement with adia- 
batic hydrodynamic simulations, but it remains unclear 
how valid it is for radiative simulations which include 
cooling, star formation, and feedback. We will take an 
alternative approach, choosing appropriate values for the 
free parameters such that they are in agreement with hy- 
drodynamic simulations, are complimentary with other 
recent semi-analytical models, and are flexible enough 
for us to apply them to our non-polytropic model. 

We set the polytropic index to be 7 = 1.2 as sug- 
gested by hydrodynamic simulations and used in other 
recent semi- analytical models (see [391 ] . and references 
therein). This value is also consistent with the range 
of values used by Komatsu and Seljak [3(1 HH, who 
parametrized the polytropic index as a function of the 
concentration parameter and found that it varies only 
weakly with c. In reality, the effective polytropic index 
7eg = dlnP g /dlnp g is likely to be a function of radius. 
Towards the center of the halo where cooling is more ef- 
ficient because of the higher densities, the temperature 
can decrease, resulting in 7 c ff < 1. In the outskirts of 
the halo beyond the virial radius where shockheating is 
less efficient, the effective polytropic index will approach 
the characteristic value 7 ff = 1.62 for the IGM 40]. Fur- 
thermore, non-gravitational feedback can also change the 
equation of state. In our second model, we consider a 
profile where the effective polytropic index is scale de- 
pendent. 



Halos found in adiabatic hydrodynamic simulations are 
generally well described by the virial theorem. The tem- 
perature at the virial radius is close to the virial tem- 
perature and the central and average temperature are 
found to be slightly higher (e.g. [JO, HH). Therefore, 
we choose to equate the temperature at the virial radius, 
T gtC = T g (c), to T v ; r and this fixes the central tempera- 
ture as, 



T, 



9,0 



3T V 



7- 



1 



7 



cf(c) 



m(c) 



(13) 



Our chosen values for T g ,o and T„ jr are also consistent 
with the range of values used by [30, HH • 

We normalize the density profile by fixing the mass 
of hot gas within the virial radius, r v ir, and compare 
the gas fractions predicted by this model to X-ray con- 
straints within rsoo. X-ray observations of hot clus- 
ters have demonstrated that the cumulative gas frac- 
tion within rsoo approaches a constant value that lies in 
the range f ffa s.Ron w 0.10 to 0.16 independent of cluster 
mass (pi| |43|). Furthermore there is good evidence that 
/ gas (< r2on ) converges to the universal baryon fraction 
fb (Hzl, HEi). These results agree with cosmological 
hydrodynamic simulations, which also suggest that for 
massive clusters, there is very little evolution of the gas 
fraction with redshift within r vir [H, [13, [H, Si, H3] . 

Observational constraints on the gas fraction in the 
cooler halos of galaxies and groups of galaxies are much 
weaker than the cluster measurements. However, there 
is evidence for a decrease of / gas (< r 50o) with decreasing 
halo mass. It also appears that / ga s(< ^200) does not 
approach the universal baryon fraction fb ([44J, see a l so 
[45j|). The main process that changes the gas fraction 
while maintaining the baryon fraction is condensation of 
cold, low entropy gas into stars. The lower value of / gas 
for smaller halos is consistent with the hig her stellar frac- 
tion measured for lower mass halos ( [5ll |52| ) combined 
with the fact that lower mass halos are more sensitive to 
non-gravitational heating which can expel gas from their 
shallower potentials. There are few observational con- 
straints on the evolution of the gas fraction with redshift 
for cooler systems, but simulations suggest that there is 
little evolution of the baryon fraction from z = 1 to z = 
down to the galactic scales of M v ; r ~ 10 12 M Q [53| . 

Taking these uncertainties into account, we assume 
that the baryon fraction within the virial radius is given 
by the cosmic fraction fif,/fi m , but allow for a fraction 
/* = 0.1 of baryons in the form of stars. We also as- 
sume that /gas^ir 

is redshift-independent (for z < 1) for 
all halo masses (M v i r = 10 13 — 10 15 M Q ) that we consider. 
The central gas density is then given by, 
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Note that Komatsu and Seljak [3Cj, [3l| chose to fix the 
gas fraction at the virial radius to the cosmic value, with- 
out any explicit allowance for stars. In general, the in- 
tegrated gas mass and pressure within the virial radius 
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FIG. 1: Integrated gas fraction within r 5 oo for the polytropic 
(solid curve) and entropy injection (dot-dashed curve) models, 
each normalized as described in the text. The model curves 
are compared to the universal baryon fraction, as measured 
from the WMAP 5 year data, with 10% removed to account 
for stars (dotted curve). The data points were obtained by 
Vikhlinin et al. [17] from a CHANDRA sample of nearby re- 
laxed galaxy clusters. 



from our model agrees well with theirs with typical dif- 
ferences of order ten percent. We note that this simple 
parametrization does indeed fit X-ray observations as can 
be seen from Fig. [T] 

Now that the density and temperature of the gas are 
completely specified for the polytropic model, we can 
write down additional quantities that are relevant to SZ 
and X-ray observations. The electron pressure and en- 
tropy are given as, 



Pe — Tl e kT e — PefiVpo\y ; 



and 



-2/3 _ c ,,7-5/3 
— J e,0i/p oly 



(15) 



(16) 



2/3 

where P e ,o = n e ^kT e fi and S e p = T e ^n e ' are the cen- 
tral electron pressure and entropy, respectively. The elec- 
tron number density n e is calculated from the gas density 
assuming a fully ionized gas and the electron temperature 
is assumed to be equal to the gas temperature. 



C. Entropy model 

Observations indicate that in the inner regions of low- 
temperature clusters there is excess entropy above the 
predictions of the self-similar model ([5J, |55j, l56j). To- 
gether with the observed departure from the simple scal- 
ing relations suggested by purely gravitational physics 
571 ] . this means that non-gravitational effects should be 
included when modeling the ICM gas distribution. Var- 
ious authors have suggested that the ICM was heated 



by some energy input, e.g. via star formation, SN ex- 
plosions or AGN feedback. For example, Voit et al. [58[ 
showed that preheating can explain the entropy profiles 
of groups. This model was also studied |59l| in the con- 
text of SZ observations. Ostriker et al. [39| constructed a 
model including various non-gravitational processes that 
could be used to match the observed scaling relations of 
clusters (see also [o0|). 

Similar to Voit et al. [HJ (see also [6l[ and [H|), we 
allow for an additive term Si n - 3 to the polytropic en- 
tropy profile that incorporates the combined effect of 
non-gravitational processes, such as feedback from su- 
pernovae and galactic nuclei, and radiative cooling and 



star formation. 



Si n j is constant with radius so that 
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The addition of an entropy term modifies the tempera- 
ture profile according to 



Tg( x ) = T gA y2nt\ x ) + T inj y\ 



2/3 



(18) 



so that the central temperature is now T g (0) = T g- i+T- m y 
Here Tj n j is the amount of injected thermal energy per 
particle, and relates to the amount of injected entropy 
via 
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keV cm 2 . (19) 



The form of heating proposed here is effective at in- 
creasing the temperature in the inner regions of the halo 
but has little effect in the outskirts of the halo for all 
but the lowest mass halos (see Fig. [T5]) . The additional 
entropy due to the non-gravitational heating term breaks 
the self-similarity of the cluster physics. This very simple 
model provides a good phenomenological description that 
is adequate for our purposes in spite of its shortcomings. 

We assume that the gas remains in hydrostatic equi- 
librium with the dark matter, which yields an implicit 
equation for the modified profile y en t, which we solve nu- 
merically. For a given mass and redshift, there are 4 free 
parameters in this model: the index 7 that character- 
izes the polytropic part of the temperature, the central 
gas temperature, the central gas density and the amount 
of injected energy Ti n j. In general the gas will not re- 
tain a polytropic equation of state and its profile will be 
altered. However, we want to continue to associate the 
polytropic-like terms with virialization and shockheating, 
thus we keep 7 = 1.2 and normalize the polytropic term 
in the temperature equation just like in the polytropic 
model. That means we set the value at the virial radius 
to be equal to the virial temperature, which fixes the 
constant T q 1 via 



T, 



9,1 



2/ent ( C ) 



(20) 



Note that we implicitly assume that the gas has first 
settled in the gravitational potential set up by the dark 
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FIG. 2: Entropy at 0.1r v i r against halo mass, M v i t /Mq, for 
the polytropic model (solid curve) and the entropy model 
(dot-dashed curve), showing a break in the scaling of the core 
entropy scaling as suggested by X-ray observations. 



matter before being redistributed by the injection of en- 
tropy. 

We normalise the gas density by assuming that the gas 
pressure at the virial radius remains unchanged by the en- 
tropy injection, i.e. P^(r y ii) = -Pg a ° ly ( r vir)- This means 
that we assume that the gas settles back into a pressure 
balanced hydrostatic equilibrium at the virial radius af- 
ter the energy injection. Note that with our choices, the 
entropy model reduces to the polytropic model when T; n j 
goes to zero. The resulting electron temperature, density, 
pressure and entropy profiles for the polytropic model 
and entropy injection model are shown in Figs. [T3]fT5l in 
the Appendix, where they are discussed in more detail. 

In order to adjust the amount of heating Tj n j, we 
construct a fitting function for T- m j such that S- m j « 
100 keVcm 2 independent of halo mass and redshift. It 
has been shown that this form and amount of feedback 
reproduces the observed scaling relations of clusters 
(see [23| for a detailed discussion of the feedback energy 
available from physical processes such as supernovae or 
AGN heating). 

In Fig. [2] we plot the entropy, S e at 0.1 r v ; r against virial 
mass for the different models. It can be seen that the 
entropy model has more central entropy in group halos. 
This results from the feedback that heats the gas and 
flattens the density profile, by pushing more gas from 
the center to the outskirts, thereby breaking the scale 
invariance of the entropy-mass relation. We also note 
that the level of feedback chosen in our entropy injection 
model provides an entropy floor of S « 100 — 200 keV cm 2 
which agrees with en trop y profiles derived from X-ray 
observations (see [H, [Hf) and simulations (see (63[) of 
galaxy clusters and groups. 

An additional constraint on the entropy injection pa- 
rameter comes from the X-ray luminosity-temperature 



FIG. 3: Integrated X-ray luminosity within r 5 oo against 
emission- weighted gas temperature, T cw , in the polytropic 
model (solid curve) and the entropy model (dot-dashed 
curve). The data points were obtained by McCarthy et al. 
[H from a compilation of CHANDRA and XMM-NEWTON 
data. 

relation. The X-ray luminosity within radius rx is 
Lx «r x) = 1 .41x 1 0« el , 8 -.( T5 ^) 2 ( j ^) 1 

x ( j^T J x 47r y yx(x p ) x p dx p 



where the projected radius is r p — x p r s — y 'r 2 — P r 2 . 
In the case of the entropy model the integrand y>x is 
given in terms of j/ gas , the dimcnsionless gas profile, as 

y x {r) = V^M + mnj/Teo)^ 3 ^) 

while in the case of the polytropic model ~$>x is simply 
obtained by setting T- m j = in the above expression. 

In Fig. [3] we plot the integrated X-ray luminosity 
within r5oo against the emission weighted electron tem- 
perature. We see that the entropy injection model 
has significantly lower X-ray luminosity for group ha- 
los, resulting in the observed break of the scale invariant 
Lx — Tx relation on group scales. 

III. MULTI-FREQUENCY FILTERING OF THE 
HALO SZ SIGNAL 

We now study the significance with which the halo 
SZ signal from nearby groups and galaxies can be de- 
tected with forthcoming multi-frequency CMB experi- 
ments. Microwave maps contain not only the thermal 
SZ, but a host of contaminants including primary CMB 
anisotropics, kinetic SZ, microwave emission from galac- 
tic dust, infrared and radio point sources, over and above 
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the detector noise of the experiment. The importance of 
utilizing the SZ as a cosmological probe has prompted 
several authors to develop specialized techniques for de- 
tecting galaxy clusters through the SZ effect. Proposed 
techniques include the maximum entropy method [&4| . 
fast independent component analysis [65| . matched fil- 
ter analysis (66|, wavelet filtering [67| and Wiener fil- 
tering [68]. Map filtering, in general, utilizes both the 
spatial and frequency information to separate galactic 
foreground emission and extragalactic point source con- 
tamination from the primary CMB and thermal SZ sig- 
nals. 

For the purposes of this investigation, which involves 
assessing the level of detection of a SZ halo, we utilize a 
simple multi-frequency Wiener filtering technique as de- 
scribed in Tegmark and Efstathiou [68[ to separate the 
halo SZ signal from other components. The method al- 
lows us to determine the level of residual noise in the 
filtered maps and a signal to noise ratio for each halo. 
Foreground and noise subtraction is done in harmonic 
space which allows one to exploit the fact that contami- 
nants such as the CMB, galactic dust and extragalactic 
point sources have power spectra that differ from that 
of the thermal SZ effect. We also include a contribu- 
tion from the thermal SZ background that has the same 
frequency dependence as the halo SZ signal. The con- 
tamination arising from the superposition of SZ sources 
along the line of sight has the potential to reduce de- 
tectability and distort observable properties of the halo 
SZ signal 

A. Multi-frequency Wiener filter 

We assume that we have microwave sky maps di (r) at 
pixel position r for M. different frequencies. The signal 
in each map originates from M components Sj (r) such as 
the primary CMB anisotropy, SZ sources, galactic fore- 
grounds and extragalactic point sources, so that 

AT(r,v) = Y / f j (v)s j (r) 

3 

where fj{y) is the frequency dependence of the jth com- 
ponent. In addition to these components each map con- 
tains detector noise rii(r) which we treat as random in 
each pixel. Then we can write our observation as 

d i (r)=F ij s j (r)+n i (r) 

where the MxAf frequency response matrix F is defined 
as Fij — J Wi(is)fj \v)dv. It is convenient to absorb the 
beam response factor into the definition of the pixel noise 
[70| which allows us to set the response coefficients w; to 
unity. 

We assume that the noise has zero mean (rii) — with 
covariance 

(ni(£)n*(£')) = (2tt) 2 5(£ - £')N l0 (l). 



We will consider the case of white noise, for which Nij{£) 
is constant. We assume that the signal and foreground 
components have means 

(si(r)) = Ai 

with covariance 

((5i(£) - A t )(S*(£') A,)) = (2tt) 2 S(£ £')S %3 {i). 

The condition that the signal and noise components are 
uncorrelated ensures that the signal covariance matrix, 
= ^ij^£,(j)i an d noise covariance matrix, Nij(£) = 
are diagonal. In what follows we will drop the 
tilde on harmonic space quantities. 

The most general linear estimator s of the signal can 
be constructed from the data d as 

s(r) = J W(r - r')d(r')dV, 

where W, is an AT x M. weight matrix. We use the flat 
sky approximation and work in harmonic coordinates. 
Some signals e.g., the primordial CMB, have a zero mean, 
Ai = 0. For a signal with non-zero mean, the condition 
of an unbiased estimator requires that 

„ M 

(3,(r)) = / Y^Wull) Sj {£) d 2 £ = A u 

3=1 

where we have used Parseval's theorem. 

The residual error in the maps from the noise and fore- 
grounds is given by 

<A, ; (r) 2 ) = ((S,(r) - s t (r)) 2 ) = J A 2 (£) d 2 £. (22) 

where 

AT M 

A?W - £ \(W t F - Ih\ 2 St 3 + rn^fN^. 

3=1 3=1 

The first term accounts for the contamination of the de- 
sired signal by other components while the second term 
measures the detector noise. A non-zero mixing arises 
when two or more components have similar frequency de- 
pendence. Requiring that the residual error is minimised 
we derive the Wiener filter weights 

Wi = S e F T [FS e F T + Nt}-\ 

To set the threshold for detection we compare the mean 
SZ signal to the residual noise and define the signal-to- 
noise ratio as 

0/N _ (Wr)) J s az (e)ede 

where the recovered signal is given by 
M 

3=1 
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FIG. 4: SZ temperature distortion profile (black curves, de- 
creasing from left to right) and integrated flux profile (blue 
curves, increasing from left to right) plotted against r/r v i r for 
the polytropic model (solid curves) and the entropy model 
(dotted curves). The panels from top to bottom show the 
SZ distortion and integrated flux profiles for halo masses 
Mvir = 10 13 ,f0 14 and 10 15 M respectively. The sharp fall- 
off that is visible in the temperature distortion profile for the 
polytropic model is a result of the Gaussian smoothing de- 
scribed in the text, and occurs at larger radius for the entropy 
injection model. 

The signal to noise ratio depends on the range of t over 
which the integration is performed. We chose l m in and 
(max to give the maximum signal to noise over this £ 
range - in practice this would correspond to applying the 
appropriate high pass and low pass filters to the recovered 
map. 

B. Halo Sunyaev-Zel'dovich Signal 

The upscattering of microwave background photons by 
hot electrons in the halo results in a projected Compton 
profile 

r m ax/fs 

Z/comp (r P ) = y c o X 2 J ysz(r)dl (23) 
o 

where the central Compton distortion is 

and the projected radius is r p — x p r s = ^Jr 2 — I 2 r 2 . In 
the case of the entropy model the projected pressure is 
given by 

y S z(r) = V^asW + (T^/Teo^W 

where y gas is the dimensionless gas profile, while in the 
polytropic model the corresponding expression is ob- 
tained by setting T U1J = 0. 
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FIG. 5: Halo spectra, s BZ (£), for the polytropic model (solid 
curves) and entropy injection model (dashed curves). The 
curves from top to bottom show the spectra for halo masses 
Mvir = I0 13 ,I0 14 and I0 15 M (blue, red and black) respec- 
tively. 

The resulting thermal SZ temperature distortion is 
given by 

AT{6) = j(x) T cmb y comp (6) (25) 

where 6 = r p /d,A, the CMB temperature is T cm b = 2.73 K 
and the frequency dependence is given by 

3(*) = X K + l} - 4 . x = GHz - ( 26 ) 

(e x — 1) 

The projected temperature distortion profiles are shown 
in Fig. [U We note that the polytropic and entropy in- 
jection models are almost indistinguishable for the most 
massive halos. In Fig. [4] we also show the cumulative SZ 
flux within radius rsz, which is given by 

rrsz/r s 

Ssz{< r S z) = Icmb g{x) X 2"7T / J/comp^p) Xp dx p 

JO 

(27) 

where 

4 x 

9( x ) = j( x ) (11^ ^ - 0.27 GJy. (28) 

While the SZ signatures of the largest clusters are nearly 
identical, it is clear from Fig. [4] that the SZ temperature 
distortions of less massive group halos are sensitive to 
feedback effects even for models which have the same SZ 
flux at the virial radius. 

The relevant quantities in harmonic space for deter- 
mining the detectability of a given SZ halo are the Bessel 
transform of the halo SZ distortion 

s S2 (£) = 2tt / ^ [AT(6)/j(x)] J Q {l6)Qde (29) 
Jo 
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where Jo is the Bessel function of order zero, and the SZ 
halo power spectrum 



2tt 



[AT(0)/j(x)} 2 J (£6)ed8 (30) 



which enters the reconstruction error if the filter weights 
are not diagonal in the SZ component i.e., (WF)i SZ ^ 

&i sz • 

The sharp cutoff in the Bessel transform, combined 
with the fact that the SZ distortion has not fallen to zero 
at # m ax; results in ringing of the profile spectrum and 
power spectrum in £-space. We therefore smooth the 
density profile using a Gaussian profile so that p gas (r) —> 
Pgas(r)e~ r /£ r max where £ is chosen to ensure that the 
total gas mass is unchanged. The smoothed profile falls 
off sharply after r max so in practice we integrate out to 
6*™° oth that is a few times larger than # max , which suffices 
to remove the ringing. The resulting spectra have slightly 
more (~ few percent) power on intermediate and large 
scales with the power going smoothly to zero at large £. 

In Fig. [5] we compare the SZ halo spectrum s sz (£) in 
the polytropic and entropy injection models for halos of 
mass M v i r = 1O 13 M at redshift z = 0.01, mass M v ; r = 
10 14 M Q at redshift z = 0.1, and mass M vlr = 10 15 M© 
at redshift z = 1.0. In the largest mass halo we note 
that both our gas models produce nearly identical spec- 
tra due to the fact that the heating is small relative to the 
thermal energy of the hot gas. For the lower mass halos 
we observe that the polytropic model produces a larger 
amplitude SZ signal due to the higher density of gas in 
the central regions of the halo. The temperature increase 
that occurs in the entropy injection models is insufficient 
to compensate for the reduced density, which results in 
a SZ profile with smaller amplitude for larger entropy 
injection. The trend is monotonic suggesting that, for 
experiments with sufficient sensitivity and frequency cov- 
erage to extract the SZ halo signal, the amount of entropy 
injection may be measurable from these SZ observations. 



1. Galactic dust emission 

We model the frequency dependence of thermal galac- 
tic dust emission as 



/dust(^) = c(x) c*(x) 



,.3+Q 

''dust 



g^-dust — \ 

where a is the emissivity index and 



£dust = hv/k B T dnst . 

(31) 



c(x) 



x 2 2k V k„T r 



he 



(32) 

convert antenna temperature to thermodynamic temper- 
ature and specific intensity to antenna temperature re- 
spectively. In our model we assume an emissivity index 
a =1.7 and a dust temperature T dust =18 K (zHzHzl- 
We model the spatial power spectrum of the thermal 
dust component as a power law 



(33) 



where (3 is the power law index. We set f3 = 3 which was 
the value derived from an analysis of the DIRBE maps 
[74| . We fix the amplitude of the galactic dust emission 
to be A = 10.2 fiK at 56.9 GHz. 

An analysis of the FIRAS and DIRBE datasets [7f| 
[76j has provided evidence for two dust components with 
different temperatures and emissivities. We account for 
uncertainty in the emissivity by introducing a residual 
dust component with the same spatial power spectrum 
but with scatter, 5a, in the emissivity index. We choose 
5a = 0.3 as suggested by the analysis of Finkbeiner et al. 
[77j . which is consistent with the results of Draine and 
Lazarian [72j. The top left and top right panels of Fig. [5] 
display the power spectra of the dust and residual dust 
components respectively. 



2. Radio and infrared point sources 



C. Foreground components 

We assume that the spatial and frequency dependence 
of each foreground component can be written as a prod- 
uct 



AT ( 2 (M) 



A 



(i)/(i)W C AB 



where /(,) (v) measures the average frequency dependence 
of a component, CV is its spatial power spectrum and Ar{\ 
Note that f{i){v) gives the frequency dependence of the 
RMS fluctuations in thermodynamic temperature refer- 
enced to the CMB blackbody. We normalise the fre- 
quency term, f{i){v), to be unity at z/» = 56.9 GHz and 
the spatial term CWj) to be unity at £* = 2 so that the 
units are absorbed into the overall amplitude Auy We 
now consider models for each of the foreground compo- 
nents in turn. 



We consider two point source populations: radio 
sources e.g., blazars, and infrared point sources e.g., early 
dusty galaxies. We model radio sources using a fit to the 
WMAP Q-band {v = 41 GHz) data [z|: 



dN 
dS v 



No 
Sn 



-2.3 



(34) 



where iV = 80 deg" and S a = 1 mJy. Since we are 
mostly concerned with fluxes at the mJy level, the slope 
of the distribution was altered to —2.3 from the fiducial 
slope of —2.7 17911 . In the case of infrared point sources 
we use the fit |80l | to SCUBA observations [8lJ at v Q = 
350 GHz, given by 



dN 
dS v 



No 

So 



3.3' 



(35) 
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FIG. 6: Angular power spectra of the foreground components 
and ACT detector noise at 145 GHz (solid curves), 215 GHz 
(dotted curves) and 280 GHz (dashed curves) respectively. 
The top left panel displays the power spectra of the galactic 
dust, thermal SZ (tSZ) background and detector noise. The 
top right panel displays the power spectra of the residual tSZ 
background and residual galactic dust. The bottom left and 
bottom right panels display the spectra of infrared and radio 
point sources respectively. All spectra are compared to the 
lensed CMB spectrum (bold solid curve) in each panel. 



TABLE I: Experimental specifications for the PLANCK |26j, 
ACT [IK] and SPT [2f| experiments. In addition to the nomi- 
nal ACT and SPT surveys, the specifications for a wider ACT 
survey and deeper SPT survey are also listed, where we have 
assumed a fixed total integration time in rescaling the pixel 
noise. 



software package 1 using the WMAP 5 year best fit cos- 
mological model. The lensed CMB angular power spec- 
trum is shown in Fig. [6] 



where N a = 1.5 xl0 4 deg" 2 and S a ~ 1.8 mJy. The 
sources fluxes were extrapolated to the frequencies of 
the various CMB experiments using a power law, S v oc 
(v/v ) a . We use a spectral index a = for radio sources 
and a — 2.5 for infrared sources {79| . 

The power spectra of the point sources is calculated 
from 

^-(i)T> (36 » 

where dB/dT is the derivative of the Planck spectrum 
and Scut is the imposed flux cut, which we assume to be 
5 mJy for ACT and SPT ^ and 250 mJy for PLANCK 
[H|]. We assume that the point sources are spatially un- 
correlated on the sky, thus the power is constant at all 
multipolcs. The power spectra of infrared and radio point 
sources are displayed in the bottom left and bottom right 
panels of Fig. [5] respectively. 



3. Cosmic microwave background 

The cosmic microwave background anisotropy has a 
constant frequency dependence with reference to the 
blackbody temperature so that f C mh{v) = 1. The CMB 
power spectrum, Ci, was calculated using the CAMB 



4- SZ background 

The projection of SZ sources of varying mass and red- 
shift along the line of sight creates a diffuse SZ back- 
ground which can contaminate halo SZ observables. To 
account for the contamination of the foreground halo 
signal by background clusters, we model the SZ back- 
ground statistically by using its power spectrum. Ide- 
ally one would utilise a simulated map of SZ halos to 
study the contamination due to projection effects but 
we defer this investigation to a future publication. This 
map would also take into account the contamination from 
hot gas outside collapsed structures, though Hernandez- 
Monteagudo et al. f83| have shown that this component 
does not significantly contribute to the thermal SZ power 
spectrum. 

The frequency dependence of the thermal SZ back- 
ground is the same as that of the thermal SZ halo sig- 
nal given in [26l The power spectrum of the SZ back- 
ground is computed following Komatsu and Seljak [3lj 
over the mass range 10 12 — 10 16 M Q . We also allow for 
an uncertainty in the SZ background which we conser- 
vatively model as arising from background halos smaller 
than 5 x 10 14 M Q . The top left and top right panels of 
Fig. [5] display the power spectra of the SZ background 



CAMB: |http: / /www.camb.info| 
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FIG. 7: The minimum detectable mass, M v i r , of SZ halos 
as a function of redshift, z, for the ACT experiment. The 
minimum detectable mass is plotted for an entropy model 
with S/N = 5 (solid curve) and S/N = 3 (dotted curve), 
and a polytropic model with S/N=5 (dashed curve) for the 
nominal ACT survey. Also shown is the minimum detectable 
mass with the ACT wide survey for an entropy model with 
S/N — 5 (dot-dashed curve). 



and residual SZ background. 



D. Detector noise and experimental specifications 

We model the detector noise as an additional sky signal 
[70| with power spectrum 



N f 



(37) 



in a given frequency band, i. In this case each of the 
sky signals is not convolved with the experimental beam. 
We assume that the experimental beam is a Gaussian of 
width 6]j t i so that the full width at half maximum is given 
by FWHM, = V81n2 6> M . The inverse noise weight 



1 = a ii x 



(38) 



is defined as the noise variance per pixel times the pixel 
area in steradians. 

We consider three nominal experiments, a shallow all 
sky SZ survey by the PLANCK surveyor 2 , a 200 deg 2 
SZ survey by the Atacama Cosmology Telescope (ACT) 
3 and a 4000 deg 2 survey by the South Pole Telescope 
(SPT) 4 . We also consider a wider ACT survey (over 4000 
deg 2 ) and a deeper SPT survey (over 200 deg 2 ) where we 



PLANCK: http:/ /www.rssd.esa.int/index.php?project=planck 



ACT: 



http: / /www. physics. princeton.edu/act 



SPT: http://pole.uchicago.edu/ 



have rescaled the pixel noise using a fixed total integra- 
tion time. Specifications for the various experiments are 
listed in Table U @, Efj. The top left panel of Fig. © 
displays the detector noise power spectra for the ACT 
experiment. 



IV. DETECTABILITY OF THE HALO SZ 
SIGNAL 



We now study the detectability of SZ halos in our mass 
and redshift range for the polytropic and entropy injec- 
tion models presented above. In Fig. [7] we plot the mini- 
mum detectable halo mass, or threshold mass, as a func- 
tion of redshift for the ACT experiment using a detec- 
tion significance of S/N — 3 and S/N = 5. Similar plots 
for the PLANCK and SPT experiments are presented in 
Fig. [5] and Fig. [5] respectively. 

At high redshift (z ~ 1) the ACT experiment reaches 
a threshold mass of 2 x 1O 14 M for a signal to noise ratio 
of five and sky coverage of 200 deg 2 , which is similar to 
the completeness limit presented in (84|. The threshold 
mass at high redshift is similar for the SPT experiment 
but larger (~ 1O 15 M ) for the all-sky PLANCK survey 
which has larger pixel noise. At low redshift, however, it 
is interesting to note that the threshold mass for ACT, 
in the case of the entropy injection model, drops below 
10 14 M Q at z fa 0.05 and as low as 4 x 10 13 M Q atz« 
0.005 (ps 20Mpc), and is similar in the case of the SPT- 
deep survey. Although the threshold mass at low redshift 
is higher for the SPT wide (4000 deg 2 ) survey and the all- 
sky PLANCK survey it is still possible to detect group- 
sized halos below 10 14 M Q at z < 0.02. This indicates 
that the thermal SZ effect in nearby group-sized halos can 
be detected with multi-frequency observations that reach 
pixel sensitivities of a few /xK, and that the detectability 
of these halos can be improved with longer integration 
times. 

For the polytropic model the detection levels are more 
optimistic because the SZ signal is larger due to the gas 
being more concentrated, however, as discussed in the 
previous section this model is less realistic, particularly 
for low mass clusters and groups. To determine the signif- 
icance with which one can distinguish the entropy model 
from the polytropic model via measurements of the SZ 
distortion in these halos we computed the \ 2 statistic 
that compares the difference in spectra between these 
models to the residual noise and foreground level from 
the nominal ACT experiment. In particular we were 
interested in the leverage one gains from the detection 
of galaxy groups and low mass clusters. In Fig. [10] we 
plot the x 2 statistic as a function of virial mass. We 
observe that the information gained from group halos 
is of the same order of magnitude as the information 
gained from larger clusters because the larger differences 
in the SZ spectra, that results from the increased impact 
of the entropy injection in lower mass halos, compen- 
sates the larger residual noise. For group halos detected 
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FIG. 8: The minimum detectable mass, M v i r , of SZ halos 
as a function of redshift, z, for the SPT experiment. The 
minimum detectable mass is plotted for an entropy model 
with S/N = 5 (solid curve) and S/N = 3 (dotted curve), 
and a polytropic model with S/N=5 (dashed curve) for the 
nominal SPT survey. Also shown is the minimum detectable 
mass with the SPT deep survey for an entropy model with 
S/N — 5 (dot-dashed curve). 
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FIG. 9: The minimum detectable mass, M v i r , of SZ halos as 
a function of redshift, z, for the PLANCK experiment. The 
minimum detectable mass is plotted for an entropy model 
with S/N = 5 (solid curve) and S/N = 3 (dotted curve), and 
a polytropic model with S/N— 5 (dashed curve). 



ing SZ cluster surveys e.g., [84[. In choosing this cutoff 
mass our aim is to quantify the additional yield of SZ 
halos, over and above the yield of more massive clusters, 
in these surveys. The anticipated number of detectable 
galaxy groups and low mass clusters are given in Table HTl 
for the different surveys. We note that the mass function 
is steep so the yield is very sensitive to the minimum and 
maximum mass limits, consequently the numbers quoted 
here should only be taken as a rough guide to the an- 
ticipated yields. We observe that all surveys will yield 
a reasonable number of detectable halos below the mass 
cutoff at the S/N = 5 level, with the numbers increasing 
significantly for halos detected at the lower significance 
of S/N — 3, though the contamination will be higher at 
this level. It is generally the case that the deep surveys 
produce a higher yield than the wide surveys for a given 
model, presumably because the mass function is so steep 
in this mass range. This trend is reversed in the case of 
the SPT wide survey at the S/N = 3 level, however, be- 
cause the increased sky area is sufficient to compensate 
for the reduced sensitivity. We also note that the yields 
for the polytropic model are larger than those for the 
entropy model due to the larger signal in the polytropic 
model. 

In Fig. [11] we compare detection curves for the ACT 
experiment to the distribution of nearby groups in the 
USGC (UZC-SRSS2 Group Catalog) [g| and a group 
catalog compiled from the Sloan Digital Sky Survey 
(SDSS) HI]. The USG catalog is based on the updated 
Zwicky Catalog (UZC) and Southern Sky Redshift Sur- 
vey (SRSS2), and contains 1168 groups of galaxies out 
to a redshift of z ~ 0.04 and over a solid angle of 4.69 
sr, while the group catalog based on the SDSS contains 
301237 groups from z = 0.02 to z = 0.2. Comparing the 
ACT detection curves to the USGC we find that there are 
429 groups and clusters in the catalog out to a redshift 
z = 0.04 that can be detected at a level S/N — 5, with 
222 of these halos having a mass below 2 x 10 14 M Q . At 
a signal to noise level of S/N = 3 the numbers increase 
to 488 halos with 281 of these halos having a mass below 
2 x 10 Mg. At redshifts beyond those probed by the 
USGC the SDSS catalog contains additional groups that 
are above the minimum detectable mass level of SZ ex- 
periments. The comparison of the threshold mass curves 
to these catalogs suggests that there are a large num- 
ber of galaxy groups already detected in redshift surveys 
that can be detected through targeted observations with 
upcoming SZ experiments. 



with higher significance these measurements will provide 
useful joint constraints on the gas fraction and level of 
entropy injection in these halos. 

It is interesting to quantify the yield of galaxy groups 
and low mass clusters that are detectable in these sur- 
veys. We calculate the yield by integrating the cluster 
abundance from the minimum survey threshold mass to 
a cutoff mass of 2 x 10 M©, which is roughly the mini- 
mum mass quoted for the detection of clusters in upcom- 



V. DISCUSSION AND CONCLUSIONS 

We have investigated the detectability of SZ groups 
using an analytic prescription for the hot gas in these 
halos. The models that we studied were based on hot 
gas being in hydrostatic equilibrium with the dark mat- 
ter halo, and described by a polytropic equation of state, 
or an equation of state modified to include an entropy 
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Experiment 


Model 


S/N = 3 


S/N = 5 


ACT Deep 


Entropy 


1520 


696 


ACT Deep 


Polytropic 


3570 


813 


ACT Wide 


Entropy 


307 


23 


ACT Wide 


Polytropic 


1371 


87 


SPT Wide 


Entropy 


5487 


131 


SPT Wide 


Polytropic 


15856 


759 


SPT Deep 


Entropy 


2323 


706 


SPT Deep 


Polytropic 


3845 


1465 


PLANCK 


Entropy 


653 


227 


PLANCK 


Polytropic 


1110 


418 



TABLE II: Anticipated yield of SZ halos detectable by the 
ACT, PLANCK and SPT experiments, calculated by inte- 
grating the cluster abundance from the minimum detectable 
mass for a given survey and S/N level to an upper mass limit 
of M vir = 2 x 10 14 M Q . 



*. -I Entropy Model - S/N = 5 
*;;Polytropic Model - S/N = 5 



10 14 - 



0.001 




1.000 



FIG. 11: The minimum detectable mass, M v -„, of SZ halos as a 
function of redshift, z, for the ACT experiment compared to a 
compilation of groups from the USGC [§5| (plotted as crosses) 
and a galaxy group catalog based on the Sloan Digital Sky 
Survey (SDSS DR4)[§i] (plotted as dots). The SDSS groups 
are only shown up to a redshift of z = 0.1. The minimum 
detectable mass is plotted for an entropy model with S/N = 5 
(solid curve) and S/N = 3 (dashed curve). 



10" 



FIG. 10: Significance with which the entropy model can be 
distinguished from the polytropic model by the nominal ACT 
survey as a function of halo mass, M v i r . The x 2 statistic is 
calculated by comparing the difference between model spec- 
tra to the residual noise and foreground level of the ACT 
experiment. 



injection term. We have found that the entropy mod- 
els are distinguishable from the polytropic models via 
measurements of their SZ distortion, even in low mass 
clusters and galaxy groups. While these models provide 
a useful starting point to evaluate the detectability of SZ 
groups, an improved analysis will include a more realis- 
tic treatment of the gas distribution as provided by high 
resolution cosmological simulations, which we intend to 
pursue in a forthcoming paper. 

Another issue that we have only partially addressed 
here, through the inclusion of an SZ background con- 
taminant, is the confusion caused by the superposition 
of SZ distortions from hot gas in halos along the line of 
sight (see for e.g., [69|, Here again a large volume 

cosmological simulation will help to quantify the impact 



of the SZ background on the detection of group halos and 
their recovered flux. By taking advantage of the fact that 
nearby group halos produce a more extended SZ signal, 
we aim to mitigate the impact of the SZ background by 
devising algorithms to separate the group signal from the 
SZ emission at smaller angular scales produced by higher 
redshift clusters. Finally the combination of maps of the 
various foreground contaminants with simulated SZ maps 
will allow us to undertake a more accurate treatment of 
the foreground contamination. While we have been rel- 
atively conservative in our modelling of the foreground 
contaminants, we have not included effects such as the 
clustering of infrared point sources, which could turn out 
to be a significant contaminant in the extraction of SZ 
halos [Hj]. We have also not included the kinetic SZ ef- 
fect as a possible contaminant because it has the same 
frequency dependence as the primary CMB but a much 
smaller amplitude on the relevant angular scales. For the 
same reason we have not attempted to detect the halo 
via its kinetic SZ signal, though we note that a detec- 
tion of the kinetic SZ signal could be enhanced by cross- 
correlation with optical or X-ray observations of the halo 
or the thermal SZ signal. 

Prospects for detection of SZ clusters have been stud- 
ied previously in the case of PLANCK gl, ACT 
[HED] and SPT [H, but these studies have mainly fo- 
cused on the statistics of SZ detections above the mass 
completeness limit of the respective surveys. Pace et 
al. [9l| found that ACT could detect SZ halos down to 
6 x lO 13 /i -1 M0 fairly independent of redshift, whereas 
we have found that these halos become hard to detect 
at higher redshifts. The analysis of [9l| only included 
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FIG. 12: Electron number density profiles, n e (10 _3 cm~ 3 ), 
plotted against r/r v i r for the polytropic model (solid curves) 
and entropy model (dot-dashed curves). The panels from top 
to bottom show the density profiles for halo masses M v i r = 
10 13 ,10 14 and 10 15 M respectively. 



FIG. 13: Electron temperature profiles, T e (keV), plotted 
against r/r v i r for the polytropic model (solid curves) and 
the entropy model (dot-dashed curves). The panels from 
top to bottom show the temperature profiles for halo masses 
M vil = 10 13 , 10 14 and 10 15 M Q respectively. 



the CMB as a foreground so it is conceivable that the in- 
clusion of point source foregrounds would degrade their 
forecasts for low mass halos at high redshift, as we have 
found to be the case in our analysis. We also note that 
the analysis of [II] found that PLANCK could detect 
halos of mass 6 x 1O 13 M //i below z — 0.1 when they in- 
cluded the CMB and all galactic foregrounds. In our 
analysis we have emphasised that smaller halos, with 
M ~ 3 - 4 x 1O 13 M , could be detected at z < 0.01. 

The detection of hot gas in these galaxy groups and 
low mass clusters with the upcoming set of SZ survey 
experiments will provide an interesting probe of galaxy 
formation and its effect on the distribution and state of 
the hot gas, which is most prominent in these halos. A 
measurement of the SZ distortion at the virial radius of 
these halos will set a joint constraint on the level of en- 
tropy injection and the baryon fraction, which can be 
compared to the predictions from galaxy formation mod- 
els. In particular a measurement of the baryon fraction 
in the outskirts of galaxy groups and low mass clusters 
will provide a unique update to the baryon census in the 
local universe. 

While the upcoming generation of SZ experiments have 
been designed to carry out blind surveys of galaxy clus- 
ters, a targeted survey of galaxy groups already detected 
in optical and X-ray observations would yield a very in- 
teresting set of objects to study. The measurement of a 
diffuse signal on the scale of tens of arcminutes will be 
challenging though, and carefully planned and executed 
observations will be necessary to control systematic ef- 
fects and produce high fidelity maps of the SZ distortion 
in these halos. In combination with existing optical and 
X-ray observations of these halos, these measurements 
will enhance our knowledge of the physics of galaxy for- 
mation and its effect on the intragroup medium. 



FIG. 14: Electron pressure profiles, P e (10 -3 keVcm -3 ), plot- 
ted against r/r v j r for the polytropic model (solid curves) and 
the entropy model (dot-dashed curves). The panels from 
top to bottom show the pressure profiles for halo masses 
Mvir = 10 13 , 10 14 and 10 15 M respectively. 
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APPENDIX: GAS PROFILES 

We present here the radial profiles of the electron tem- 
perature, density, pressure and entropy for the polytropic 
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FIG. 15: Entropy profiles, S e (keVcm 2 ), plotted against 
r/rvir for the polytropic model (solid curves) and the entropy 
model (dot-dashed curves). The panels from top to bottom 
show the entropy profiles for halo masses M vir = 10 13 ,10 14 
and 10 15 Mq respectively. 




FIG. 16: Integrated X-ray luminosity, Lx(< r/rvir), plot- 
ted against r/r v i r for the polytropic model (solid curves) and 
the entropy model (dot-dashed curves). The panels from top 
to bottom show the integrated luminosity profiles for halo 
masses M vil = 10 13 , 10 14 and 10 15 M respectively. 



model and entropy injection model (see Figs. fT3]fT5|) . For 
the higher mass halos (~ 1O 15 M ), the entropy injec- 
tion model profiles are very similar to the corresponding 
polytropic model profiles, which indicates that the distri- 
bution of gas in large clusters is fairly insensitive to the 
injection of entropy. 

For lower mass halos, the imposed heating is much 
more effective, resulting in a higher electron temperature 
especially in the central parts of the halo (see Fig. [T3|) . 
This reflects the fact that feedback effects are more sig- 
nificant for galaxy and group sized halos, raising the tem- 
perature above the shock-heated infall value. The density 
profiles of the polytropic and entropy injection models are 
significantly different for the low halo mass range, demon- 
strating that the entropy injection has a more marked 
effect on group sized halos, pushing gas into the outer re- 
gions of the halo and flattening the density profile. There 
is much more hot gas in the inner halo regions in the poly- 
tropic model which produces levels of X-ray emission in 
galaxy sized halos that are in violation of observational 
constraints, as we discuss below. Similarly, the imposed 
heating only significantly alters the pressure profiles for 
galaxy and group sized halos (see Fig. Q3]), greatly lower- 
ing the central electron pressure. The heating term was 
modelled such that the pressure was unaltered at the 
virial radius, where we do not expect feedback to have 
an effect even for the lowest mass halos. The entropy in- 
jection model was constructed to have significantly more 
entropy than the polytropic model in the inner regions 
for the low mass halos, where the injected energy input 
is significant relative to the gravitational binding energy 
of the halo. 

The cumulative X-ray luminosity profiles for the two 
models are shown in Fig. [T|)J While the integrated lumi- 
nosity profiles are similar for high mass halos, the X-ray 
luminosity in the polytropic model is significantly higher 
for low mass halos, in violation of the upper limits on the 
diffuse X-ray emission from hot gas in nearby galaxy ha- 
los like M31 @,[92]. The entropy model does not violate 
these constraints though, as the reduced central density 
in this model lowers the X-ray luminosity in the inner 
regions despite the increased temperature. 
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